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ABSTRACT 

Aims. We discuss instrumental and analytic methods that have been developed for the first generation of bolometric cosmic microwave 
, background (cmb) polarimeters. The design, characterization, and analysis of data obtained using Polarization Sensitive Bolometers (PSBs) 
' are described in detail. This is followed by a brief study of the effect of various polarization modulation techniques on the recovery of sky 

polarization from scanning polarimeter data. 

Methods. Having been successfully implemented on the sub-orbital Boomerang experiment, PSBs are currently operational in two terrestrial 
CMB polarization experiments (Quao and the Robinson Telescope). We investigate two approaches to the analysis of data from these 
experiments, using realistic simulations of time ordered data to illustrate the impact of instrumental effects on the fidelity of the recovered 
polarization signal. 

Results. We find that the analysis of difference time streams takes full advantage of the high degree of common mode rejection afforded by 
the PSB design. In addition to the observational efforts currently underway, this discussion is directly applicable to the PSBs that constitute the 
polarized capability of the Planck HFI instrument. 

I 
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1. Introduction 

Recent advances in millimeter-wave instrumentation and techniques have transformed observational Cosmic Microwave 
Background (cmb) research. Over the course of the past decade, statistical detections of the minute temperature variations in 
the CMB have given way to high signal to noise imaging of the surface of last scattering. 

The rapid pace of technological development in the field of cmb research has contributed to an equally remarkable rate of 
progress in our understanding of the Universe; cosmology is in the midst of an abrupt transition from a data-starved theoretical 
framework to a rigorously tested standard model. Remarkably, as of early 2006, the currently available cmb data are sufficiently 
precise (and accurate!) that, within the framework of the most simple inflationary models, the majority of the scientific potential 
of the cmb temperature anisotropics has already been realized. Measurements of the polarization of the cmb provide an important 
confirmation of the validity of the theoretical framework through which the data are interpreted, and also aid in the precise 
determination of the parameters of the theory. 

Further motivation comes from the tantalizing possibility that the cmb polarization holds a unique imprint from gravitational 
waves generated during the epoch of Inflation. A detection of this signature would represent a probe of physics beyond the 
standard model — both that of cosmology, and potentially that of particle physics. 

Current experimental efforts are focused on the development of high fidelity cmb polarimeters capable of characterizing 
the small fractional polarization of the cmb. The Boomerang03 experiment was the first bolometric instrument to measure the 
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Fig. 1. The corrugation geometry of the Boomerang PSB feeds, with dimensions in millimeters. The radiating aperture is on the right hand 
side, while the PSB module is seated on the left. 

polarization in the cmb, and the first of several to use the Polarization Sensitive Bolometers (PSBs) developed for Planck HFI. 
Two additional terrestrial telescopes using PSBs, Quad and the Robinson Telescope, are currently observing from the South Pole. 

We describe the experimental approach employed by the first generation of bolometric instruments to successfully probe cmb 
polarization and discuss aspects of the data analysis which may inform future observations. Following a brief description of PSBs, 
we provide a pedagogical description of the method of analysis that has been applied to the B03 data ll27ll34ll35Ll3q. 14311 . Finally, 
we investigate the merits of several modulation schemes for scanning polarimeters which are directly applicable to current and 
proposed cmb polarization experiments at the South Pole and from Antarctic Long Duration Balloon flights. 

2. Polarization Sensitive Bolometers 

Quasi-total power and correlation receivers (both heterodyne 0] and homodyne jl^l) using low noise front-end amplifier blocks 
based on HEMT amplifiers are mature technologies at millimeter wavelengths. The fundamental design principles of these 
receivers are well established and have been used to construct polarimeters at radio to mm-wave frequencies for many years J3 
I32ll40ll46[l4^ . Although cryogenic bolometric receivers achieve much higher instantaneous sensitivities over wider bandwidths 
than their coherent analogs, the intrinsic polarization sensitivity of coherent systems has made them the technology of choice for 
the first generation of CMB polarization experiments. 

Prior to the release of the Boomerang03 results IztI |^3||. all^ published detections of cmb polarization were obtained 
from experiments relying on the proven HEMT technology f3V20', 29, 31, 42, "4^. While interferometry is a robust method of 
polarimetry, the A^^ scaling of the complexity of correlators prohibits scaling of the design to large optical throughput, limiting the 
raw sensitivity of a practical interferometric experiment. Although much progress has been made in increasing the sensitivity, and 
decreasing the footprint, of single-dish HEMT based correlation receivers 1 13], this technology has not yet been demonstrated 
with the sensitivity or scalability of contemporary low-background bolometer arrays at frequencies above ~ 90 GHz. 

We briefly describe a bolometric system that combines the sensitivity, stability, and scalabihty of a cryogenic bolometer with 
the intrinsic polarization capability traditionally associated with coherent systems | 28]. In addition, the design obviates the need 
for orthogonal mode transducers (OMTs), hybrid tee networks, waveguide plumbing, or quasi-optical beam splitters whose size 
and weight make fabrication of large format arrays impractical. Polarization sensitive bolometers (PSBs) are fabricated using 
the proven photolithographic techniques used to produce "spider web" bolometers, and enjoy the same benefits of reduced heat 
capacity, negligible cross section to cosmic rays, and structural rigidity |52]. Finally, unlike OMTs or other waveguide devices, 
these systems can be relatively easily scaled to ~ 600 GHz, limited at high frequencies only by the ability to reliably manufacture 
sufficiently small single-moded corrugated structures. Receivers using this design have been demonstrated at 100, 150, 217 and 
353 GHz. 

Polarization sensitivity is achieved by controlling the vector surface current distribution on the absorber, and thus the effi- 
ciency of the ohmic dissipation of incident Poynting flux. This approach requires that the optics, filtering, and coupling structure 
preserve the sense of polarization of the incident radiation with high fidelity. A multi-stage corrugated feed structure and coupling 
cavity has been designed that achieves polarization sensitivity over a 33% bandwidth. 

The PSB design has been driven by the desire to minimize systematic contributions to the polarized signal. Both senses 
of linearly polarized radiation propagate through a single optical path and filter stack prior to detection, thereby assuring both 
detectors have identical spectral passbands and closely matched quantum efficiencies. 
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Fig. 2. A photograph of a 145 GHz Boomerang PSB absorber The diameter of the grid is 2.6 mm, while the absorber leg spacing, g, is 108 
jim. Each leg is 3//m wide. This device is sensitive to incident radiation polarized in the vertical direction due to the metalization of the Si3N4 
mesh in that direction. The horizontal Si3N4 beams evident in the photo are not metalized, and provide structural support for the device. The 
thermal conductivity between the absorber and the heat sink is dominated by the metallic leads running to the thermistor chip. 

Two orthogonal free-standing lossy grids, separated by ~ 60yL(m and both thermally and electrically isolated, are impedance- 
matched to terminate a corrugated waveguide structure. The physical proximity of the two detectors assures that both devices 
operate in identical RF and thermal environments. A printed circuit board attached to the module accommodates load resistors 
and RF filtering on the leads entering the bolometer cavity. For Boomerang, Planck HFI, QUaD, and the Robinson Telescope, 
the post-detection electronics consist of a highly stable AC readout with a system 1 // knee below 30 mHz lisl|30l|50|l. Unlike 
coherent systems, this low frequency stability is attained without phase switching the RF signal. 

We have designed the optical elements, including the feed antenna and detector assembly, to preserve sky polarization and 
minimize instrumental polarization of unpolarized light. To this end, the detector has been designed as an integral part of the 
optical feed structure. Corrugated feeds couple radiation from the telescope to the detector assembly. Corrugated horns are the 
favored feed element for high performance polarized reflector systems due to their superior beam symmetry, large bandwidth, 
and low sidelobe levels. In addition, cylindrical corrugated feeds and waveguides preserve the orientation of polarized fields with 
higher fidelity than do their smooth- walled counterparts. 

The coupling structure, which is cooled to below 1 K, ' consists of a profiled corrugated horn, a modal filter, and an impedance- 
matching section that allows efficient coupling to the polarization sensitive bolometer (see Figure^. In addition to a reduction in 
the physical length of the structure, the profiled horn provides a nearly uniform phase front that couples well to the other filters 
and optical elements in the system. The modal filter isolates the detectors from any unwanted higher order modes that may be 
excited in the thermal break. In addition, this filter completely separates the design of the bolometer cavity from that of the feed, 
which couples to the optics. The impedance-matching section (the re-expansion at the left side of Figure^ produces a uniform 
vector field distribution with a well-defined guide wavelength^ and characteristic impedance over a large (~ 33%) bandwidth. 

The detector assembly is a corrugated waveguide that is operated well above cutoff. The bolometers, which act as an 
impedance-matched termination of the waveguide cavity, are coupled via a weak thermal link to the temperature bath. The 
electric field in the cavity drives currents on the surface of the absorber, resulting in ohmic power dissipation in the bolometer. 
This power is detected as a temperature rise measured by means of matched Neutron Transmutation Doped Germanium (NTDGe) 
thermistors |6]. The bolometers each couple to a single (mutually orthogonal) linear polarization by precisely matching the ab- 
sorber geometry to the vector field of the coupling structure. The coupling structure has been tailored to ensure that the field 
distribution resulting from a polarized source is highly linear at the location of the bolometer 

Because the absorber geometry influences the field distribution within the coupling structure, a treatment of the bolometer 
cavity as a black-body is in general not valid. An important consequence of this fact is that any attempt to model an analogous 
multimoded optical system must consider interference terms between modes when calculating coupling efficiencies or simply 

' The PSB feeds in Boomerang03, Quad, and the Robinson Telescope operate at 240 mK, while the Planck HFI focal plane is cooled to 
~ 100 mK. 

^ The guide wavelength, A^, is typically 20% larger than free space, and dlog(/lg)/dlog(v) remains small over the entire range of operation. 
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trying to predict radiation patterns. The amplitude and phase of any higher order modes capable of propagating to the bolometer 
depend on the details of both the excitation and structure. Therefore, any numerical calculation would be susceptible to a large 
number of uncertainties associated with the appropriate boundary conditions at the bolometer. For this reason, it may prove 
difficult to extend the general single mode PSB design to a multimoded application without sacrificing crosspolar performance. 



3. Analysis 

3.1. Polarization formalisms 

The two most commonly used conventions for treating polarized radiation are the Jones and the Stokes/Mueller formalisms. The 
primary difference between the two approaches is that the Stokes/Mueller formalism manipulates irradiances, and therefore is 
applicable only to incoherent radiation. On the other hand, the Jones formalism models optical elements with matrix operations 
on the (complex) field amplitudes, making it the appropriate approach for coherent analysis. While the Jones formalism is rather 
intuitive, the Stokes formalism is more naturally suited to cmb analysis. In the following we introduce both approaches at an 
elementary level, and describe the correspondence between the two. A more detailed description of each approach may be found 
in Hecht U^, Jones ElElEl, Mueller j|37|l and Hamaker & Bregman fll . 

The general action of linear optical elements can be described in terms of the relationship between the input and output 
electric field vectors. The Jones matrix of an optical element is defined in terms of its action on the incident fields, 

e/ = Je,- , 

where the Jones matrix, J, of the system is a general product of the matrices describing individual components in the system. 

(1) 
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Of course, all such components may be rotated with respect to one another with the usual rotation matrices, 

J'=RJR^, 

with 

^ _ / cos iff - sin i/r 
~ y sin 1^ cos (ft 

This formalism allows a fairly complicated optical system to be described by a single matrix, which need only be derived once 
from the constituent components. 

As an example we describe an imperfect polarizer oriented at an angle if/ with respect to the basis in which the fields are 
defined. Such an object may be represented by the Jones matrix: 



J] 
6 



7] COS i// + 6 sin tff {?] - 6) cos iff sin iff 
(rj - 6) cos if/ sin if/ rj sin^ if/ + 6 cos^ if/ 



(2) 
(3) 



where rj > 6. A perfect polarizer would have rj = I and 6-0. After generous application of trigonometric identities, one recovers 
the general Jones matrix for an imperfect polarizer oriented at an angle if/ 



{rj + 6) + (t] - 6) cos 2if/ (rj - 6) sin 2if/ 

{rj - 6) sin 2ifj {rj + 6) - {rj- S) cos 2ifj 



(4) 



Each detector in a polarization sensitive bolometer pair acts as just such a partial polarizer, followed by a total power detector. 
The Stokes parameters are defined in terms of the electric field as follows: 



/ 

Q 
U 
V 



{E,E* + e,e;) 

{E,E* - EyE;} 

{e,e; + e,e;) 

i{E,,E*y - EyE*} 



<|£,|2) + <|£,,|2) 
<|£,|2) - <|£,,|2) 

2 (l^v-Eyl COS((^., - (f>y)) 

2 {\E,Ey\ sin(0, - 0,)) 



where the brackets, ( ), represent a time average and the fields are specified in a coordinate system fixed with respect to the 
instrument. For Thomson scattering of electrons in a quadrupolar radiation field there is no mechanism for the introduction of a 
relative phase between the two polarizations. Therefore, the cosmological Stokes V parameter is presumed to be zero. 
The action of linear optical elements on a Stokes vector, s, can be described in terms of the elements' Mueller matrix. 



Sf = Msi. 
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Given the definition of the Stokes parameters, one can derive the relationship between a Jones matrix, and the corresponding 
Mueller matrix. Following Born & Wolf 01 we find 

M,,. = itr((riJ(rjJ"'") , (5) 
where the cr, are the Pauli matrices: 

ij ^Q-^^ = [o -ij 

(6) 

1 \ / -/ ' 



cruder, =\ ^ cry = cr2 = y. ^ 

Applying a moderate amount of algebra to Equations |5] and |3 we find the first row of the Mueller matrix Mp for a partial 
polarizer This defines the total power detected as a function of the incident /, Q, U, and V parameters: 

Mn = + (7) 

MiQ = ^(ri^-S^) cos2(A (8) 

Mw = ^(7]^-S^) sinli/, (9) 
M/y = . (10) 

The signal from a total power detector is proportional to the Stokes / parameter of the incident radiation. Modeling a polar- 
ization sensitive bolometer as a partial polarizer followed by a total power detector, we find (ignoring, for the moment, the effects 
of finite beam size and frequency passband) the data may be expressed as a sum 

di = |f(l + e) ■ / + (1 - e) ■ (ecos2i/r,- + f/ sin 2.A/)] + «; , (H) 

where we have defined the polarization leakage term, e, such that (1 - e) is the polarization efficiency^, ifr is the orientation of 
the axis of sensitivity of the PSB, and s is the voltage responsivity of the detector For Boomerang, the value of the crosspolar 
leakage is typically ~ 5%, and less than 3% for second generation devices designed for the Planck HFl, BICEP, and QUAD. 

It should be noted that the noise contribution, n, is overly simplified in Equation ^2 See Appendix IaI for a more detailed 
discussion of the noise properties of bolometric receivers, which explain the general features of the noise power spectra shown 
in Figures |S]and|5] 



3.2. Polarized beams 

The angular response of an instrument can be characterized by the copolar and crosspolar power response functions P\\{r, 9, <p) 
and P±{r, 6, <p). In the time reversed sense these can be thought of as the normalized power at any point in space resulting from a 
Unearly polarized excitation produced by the feed element in the focal plane. That is, for a given polarization p - (||, ±), 

Pn{r,e,(h) = ^ ' . (12) 

^ |£||(r,0,0)|2 ^ ' 

For a single moded system, Pp has nothing to do with the properties of the detector Due to the presence of the modal filter 
in the throat of the coupling feed (see Figure^, the beam is a function only of the feed geometry and the optical elements of the 
system. To fully characterize the system, the polarized beam patterns must be considered separately from the detector. 

The exact definitions used for the polarizations on the sphere p - (||, ±) vary in the literature, but the standard is Ludwig's 
Third definition |33|. In any case, for small angles away from the beam centroid, they are very nearly equal to the Cartesian 
definition. 

Qualitatively, P\\ is similar to an Airy pattern; a Gaussian near the beam centroid, with a series of side-lobes. Pj^ is zero on- 
axis for most optical systems and, for on-axis systems, is also zero along the E- and H-planes. The peak of the crosspolar pattern 
is typically at the half power point of the copolar beam, and the peak amplitude relative to the copolar beam is fundamentally 
related to the asymmetry of the copolar beam. For an azimuthally symmetric system, such as a feedhorn antenna, this produces 
lobes in the 45° plane in all four quadrants of the beam. For an off-axis reflector such as the Boomerang telescope, the azimuthal 



' That is, in terms of the elements of the Jones matrix for an imperfect polarizer, the leakage e = 6^ Irf. This is the ratio between the minimum 
and peak power response to a pure linearly polarized source, which is a directly observable property of the PSB. 
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symmetry is lost and the lobes are bimodal. It is worth noting that the polarized beams generally depend on frequency as well as 
the field distance. 

A PSB detects the convolution of the polarized sky with the polarized beam, integrated over the frequency bandpass, and 
subject to the polarization efficiency of the detector.'* In the flat sky approximation, a time domain sample of a single detector 
within a PSB pair, di, may therefore be written as the sum of a signal component 



(13) 



di^^J flfv^F, jj dQ. (P||(r,) + P^(r,)) / + 7 nh) [q cos 2iA,- + U sin 2iA/ ) 

and a noise contribution. Here, the Stokes parameters are defined on the full sky and the integration variable is r, - hi - r, for 
a vector, n,, describing the pointing at a time sample, /. We have also defined the beam solid angle Q-b - ff dQ. (P\\ + P±). The 
normaUzed beam response and the polarization efliciency are given by, 

= ^ r-lrf- (14) 

For B03, the angle tp is modulated by sky rotation and the motion of the gondola. The calibration factor, s, converts the brightness 
fluctuations in /, Q, and [/ to a signal voltage. 

By rearranging Equation and dropping both the explicit spatial and frequency dependencies, the relation can be written 
more intuitively, 

*4/*,.= f.//</a[,.,P(ec„s2,,.f,si„2«J. ,15, 

We have made the simplifying assumption that we may remove the beam and polarization efficiencies from the integral over the 
sky, and then absorb these prefactors into a redefinition of the calibration constant, s - s' J dv (1 + e). 



3.3. Signal and Noise Estimation 

A great deal of effort has been devoted to the development of algorithms designed to estimate the signal and noise from noise- 
dominated data, and a rich literature has developed around the topic (for some recent examples, see Amblard & Hamilton 
f3\, Dore et al. fviS, Jarosik et al. |22]). In this section we outline in pedagogical detail the method used to estimate the signal 
and noise from the published Boomerang03 data 

An estimate of the instrumental noise properties that is both precise and accurate is required in order to avoid the introduction 
of a bias to the estimate of the power spectrum of the signal. For B03, the high signal to noise ratio of the data in fact complicates 
the noise estimation procedure. We solve for the noise and signal simultane ously using an iterative procedure adapted from that 
applied in the analysis of the data from the 1998 flight of Boomerang | 38, 45, 47]. The B03 data are unique among contemporary 
CMB experiments in that not only are the temperature maps signal dominated at angular scales approaching the beamsize, but 
the time ordered data are also characterized by signal to noise ratios of order unity. B03 is the only polarized dataset that is 
comparable in this regard to that anticipated from the Planck HFI. 

Assuming that the data, d, are well described as the sum of a sky signal and a noise contribution, d - Am + n, where A^ is 
the pointing matrix which maps time domain samples to pixels on the sky. If the statistical properties of the noise contribution 
are piecewise stationary with a (circulant) noise covariance matrix, N, defined as 



N,f = 2 («,+/ - <n)) («,+,' - <n)) , 



then the least squares estimate of the map is given by, 

= (A^N^'a)"' A^N^'d. (16) 



m 



The iterative procedure begins with the assumption of a white noise power spectrum {i.e., diagonal N), in which case Equation 
1 161 corresponds to a simple average of the data falhng in a given pixel. The noise contribution used in a given iteration on the 
solution to Equation^]is obtained from the estimate of the signal obtained in the previous iteration, 

= d - A (17) 



* This treatment is actually more general; it holds for any receiver that can be characterized as a total power detector preceded by an imperfect 
polarizer. That is, any receiver that is well described by a Jones matrix of the type 
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In the discussion that follows, keep in mind that e = ^. 
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Fig. 3. The fractional change in the cut-sky angular power spectra (pseudo-Cf) for temperature (torn panel), and polarization (bottom panels) 
from iteration to iteration of the mapmaker. The largest scales take the longest to converge, and polarization takes longer to converge than does 
temperature. Note that the scale of the polarization axes are increased by a factor of two relative to the temperature. 



Piecewise stationarity of the noise allows subsequent convolutions of d to be performed in the Fourier domain. 

For most terrestrial telescopes, which suffer from relatively high backgrounds and large atmospheric signals, the time stream 
is dominated by noise. Therefore, a good approximation to the noise covariance matrix appearing in Equation[^is provided by 
the power spectrum of the raw data. For orbital and balloon-based cmb experiments like Boomerang, the time ordered data are 
not noise dominated, greatly complicating an accurate determination of the noise. In general, an in-situ estimation of the noise is 
required due to the influence of atmospheric emission, unpredictable backgrounds, and scan-synchronous effects. As a result, the 
simultaneous estimation of both the signal and noise is required. 

An iterative solution for both N and in is possible by using an adaptation of the Jacobi method. The Jacobi method is an 
iterative approach to the solution of a general linear system of equations, such as Equation[^ that does not require the inversion 
of large matrices. The application to the solution of Equation[^is derived in Appendix IbI This algorithm is naturally suited to 
the problem of noise estimation, as the signal subtraction is an integral part of the iterative procedure. 

By iterating on the noise covariance matrix, N, as well as the signal, in, one approaches a general least squares solu- 
tion for both. This procedure has been used in the noise estimation of previous experiments that probed the cmb temperature 
anisotropics |45, 49|. In this application, the approach has been extended to a polarized data set. 

As described in AppendixiBl each subsequent iteration on the solution to Equation[^ m^+i, is calculated from the previous 
solution according to the procedure 

where 

= a ■ diag(AX'Ar'AX'(d-Affi<t), 

and the relaxation parameter, a < 1, is tuned to optimize the speed of convergence. Recall that, in the case of polarized data, 
the quantity d represents the left hand side of Equation|^ Therefore, the calculation of the matrix diag(A^N^ ' A)"' involves the 
inversion of the polarization decorrelation matrix on the right hand side of Equation|2] The great advantage of this method is 
that the convolution of the data with the inverse noise correlation matrix. 



n^+i = N^'(d- Ami) , 



(18) 
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Fig. 4. The convergence criteria for the iterative procedure is determined by a threshold on the rms amplitude of the correction. Histograms of 
the corrections to the I (left panel) and Q (right panel) maps (the change in pixel value between subsequent iterations) are shown for 5, 10, 20, 
40 and 80 iterations. 



can be efficiently calculated in the Fourier domain, without needing to invert the full time domain correlation matrix, N. This 
operation is simply the application of a Fourier filter to the signal-subtracted time stream using the inverse of the noise power 
spectrum as the filter kernel. 

The unbiased estimation of power spectra relies crucially on the ability to accurately model the noise properties of the 
instrument |8, 21]. In order to treat the noise in a self-consistent fashion as the realization of a Gaussian random process, it is 
necessary to measure and store the observed auto- and cross-correlations for all channel permutations, and each noise stationary 
subset. In the North American analysis of B03,^ the data are divided into 215 noise-stationary subsets, refeiTed to as chunks, each 
of which consist of approximately one hour of data. For each of these chunks the 36 (complex) auto- and cross-power spectra 
are calculated, binned logarithmically, and stored to disk. When used to generate noise realizations or construct filtering kernels, 
these binned spectra are interpolated to the discrete frequencies required by each subset of the data. 

The Boomerang readout electronics are AC coupled at ~ 6 mHz, and therefore there is no useful information in the time 
stream on timescales longer than that set by the stationarity of the noise. Dividing the data into these hour long subsets represents 
a tradeoff between sample variance and stationarity in the accuracy of the noise estimate. The non-stationarity of the B03 noise 
is illustrated in Figure |5l 

The chunk boundaries are chosen to maximize the accuracy of the noise estimate. The length of these chunks introduces 
a practical limit to the length, A^^, of the kernel applied in Equation^] The computational scaling is thus NdlogiNr) for each 
iteration of the mapmaker The memory requirement is also set by the degree of noise stationarity; the algorithm only requires 
the pointing and bolometer data for an individual chunk to be held in memory at any given time. For B03, the contributions of 
file writing and Fourier transforms to the run time are approximately equal, depending on the avaliable memory.'' 

The Fourier approach to the analysis requires the data within a chunk to be continuous and well characterized by a given noise 
power spectrum. About 7% of the B03 time stream is contaminated by transient events (primarily cosmic ray hits and calibration 
lamp pulses). These gaps are flagged, and replaced with fake data that are statistically consistent with the remainder of the chunk. 
The signal subtracted data are easily filled with any reasonable realization of the noise. Due to the small fraction of the data 
which are contaminated, the exact method of gap-filling has negligible impact on the final signal and noise estimates. 

For each chunk all [Nch(Nch - l)/2 -i- Nch] auto and cross power spectra are derived from the signal-subtracted time stream, 
n, obtained from Equation^]using the maximum likelihood maps derived from the full set of data. The noise spectra obtained 
in this manner are generally biased due to the effect of pixelizing the (continuous) sky signal, as well as the finite signal to noise 
with which the sky signal, in, is recovered 101 • 

Given a sufficiently high resolution pixelization, the signal variation within a pixel can be made to be negligibly small 
compared to the noise in the map. We pixelize the sky using the HEALPix method, at a resolution which coiTesponds to a pixel 
size of ^ 3.4' 1 12]. At this resolution we find the effect of pixelization to be well below the instrumental noise per pixel of the 
B03 data. As described in Masi et al. llssll and i27ll . the B03 cmb data are divided into a shallow and deep field, the latter being a 
subset of the former The noise per pixel of the deep field is roughly three times lower than the shallow field. 

The impact of the noise in the signal estimate is found to be significant for the data that constitute the shallow region of the 
B03 target field. The raw sensitivity of the instrument ultimately determines the signal to noise ratio of the time ordered data 

' The Boomerang team implemented two independent analysis of the time ordered data from the 2003 Antarctic Long Duration Balloon 
flight. The results from both pipelines are reported in Jones et al. |27],Masi et al. 1 35], Montroy et al. |36], Piacentini et al. |43]. 

* The Jacobi solver implemented by the North American Boomerang team requires ~ 120 MB of RAM and produces a converged GLS 
estimate of the signal and noise at a rate of 10 processor- i/channel/hour of data sampled at 60 Hz, when running on a 2 GHz AMD Athlon64 
X2 workstation. 
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Fig. 5. The time dependence of the (signal subtracted) noise power spectra of the Boomerang science channels as determined from the in-flight 
data. Each frame shows the power spectrum of each noise stationary subset (chunk) from a particular channel. The series of lines above 10 
mHz corresponds to the harmonics of the scan frequency. The signal band extends from 0.05 to 5 Hz. The diurnal dependence of the 1// knee 
is evident. The B345Z channel exhibited noise whose properties were neither stationary nor Gaussian, which is manifest in the low frequency 
contribution. 



and, when combined with the distribution of integration time on the sky, the fidelity of the recovered signal estimate, in. The 
error in the signal estimate, in, introduces a bias to the estimate of the noise power spectrum 0]. This bias is generally frequency 
dependent because of the finite bandwidth of the signal. The bias in the noise estimation varies from chunk to chunk as a result 
of the variation in signal to noise ratio in different parts of the map. 

The origin of this bias can be understood through closer examination of the signal-subtracted time stream, n, that is obtained 
from the estimate of the Stokes parameter maps, in, namely, n = d - Am. The data are assumed to consist of the sum of a pure 
signal and noise, d = s -i- n, giving 

n = s + n - Afn 

= n-n (19) 

where we have defined the projection of the signal error to the time stream as n s A(in - m). The raw noise power spectrum, 
(im' ), which is estimated from Equation^]differs from the true noise power spectrum, (nn' ), by the factor 

(20) 



<nnt) <nn-|-) 

For B03, the projection of the map errors to the time domain is highly correlated with the true time domain noise, and therefore the 
cross-correlation term dominates in Equation|20| The raw noise power spectra therefore tend to underestimate the true amplitude 
of the noise at frequencies within the signal bandwidth. The amplitude of the bias term in Equation|5n]is as high as 10% for the 
most poorly covered regions in the shallow field, and is below 1 % for the 175 chunks of the deep field. 

To correct for the bias present in the B03 noise estimates, we generate an ensemble of signal and noise simulations using 
a fiducial noise power spectrum and run the noise estimation procedure on each realization. The transfer function of the noise 
estimation procedure is then obtained by comparing the ensemble average of the estimated noise power spectra to the input 
power spectra. The size of the ensemble is determined by the required reduction of the sample variance at the lowest frequencies 
of interest; we find that the transfer function is characterized at the sub-percent level with seventy-five realizations. This bias 
transfer function is then used to correct the spectra obtained for each chunk of the time ordered data. A comparison of bias 
transfer functions that are typical of data in the high and low signal to noise regimes is shown in Figure|6l 

In order to produce noise realizations which accurately reflect the statistical properties of the instrumental noise, we require a 
framework in which to treat noise correlations between detectors in the time domain. Noise correlations in the data are expected 
both from fundamental considerations (see AppendixIXIi. as well as from the presence of correlated thermal/optical fluctuations, 
and crosstalk in the readout electronics. The measured noise from a given channel, 7i*, is modeled as the sum of an intrinsic 
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Fig. 6. Left panel: the power spectral density, in cmb units, of the noise for a representative chunk of the deconvolved B03 time ordered data. 
The black line is derived from the raw (signal plus noise) data, whereas the red line is the estimate of the signal-subtracted PSD. The scan 
frequency for this chunk appears at 12 mHz, and the cmb dipole (which appears as a triangle wave at the scan frequency) has been subtracted 
from the TOD prior to the noise estimation. Right panel: the amplitude of the noise bias as determined from an ensemble of signal plus noise 
simulations. The blue line is representative of the bias in a typical high signal-to-noise chunk, whereas the red line is the most extreme example 
found in the low signal-to-noise regime. Further discussion of the noise in bolometric detectors can be found in Appendix IaI 



(uncorrelated) component, «*, and the contributions from the intrinsic noise of the other channels, filtered through a (frequency 
dependent) crosstalk transfer function 

where, by definition, the intrinsic noise at each frequency is distributed as an uncorrelated Gaussian distribution. 
The observable quantities 

are the {NchiNch - l)/2 + Nch) auto- and cross-correlations of the signal-subtracted time streams, which are estimated directly 
from the time ordered data. After correcting for bias, the P/k are used to generate reaUzations of the noise time streams which 
exhibit the same correlation structure observed in the data. These noise realizations are constructed bin-by-bin in the Fourier 
domain. For each discrete frequency, we calculate the Cholesky factorization, H(/), of the complex (Hermitian positive definite) 
channel correlation matrix, 

P(/) = H(/)H^(/). 

Independent realizations of white noise are generated for each channel. Simulated data with the proper correlation structure 
are obtained by operating on the transform of these realizations with the Nd, x Nd, matrix Hikif) for each frequency bin in a 
given noise-stationary subset of the data. Once all of the frequency components are calculated, the inverse transform provides a 
correlated noise time stream for each channel that is used in the Monte Carlo pipeUne. 



3.4. Polarized mapmaking 

Estimates of the /, Q, and U parameters can be recovered by generating orthogonal linear combinations of the data. For each 
sample, i, of a given detector and a measurement of the projection of the orientation of that detector on the sky, i^,, one can 
construct the polarization decorrelation matrix defined by. 



' di ' 


f 


dijci 




. dijsi , 
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1 

yci 



yci 

2 2 



ysi 

2 2 









Q 







(21) 



where y = is a parameterization of the polarization efficiency. For simplicity we have abbreviated the trigonometric func- 
tions, whose argument is lifjj. 

In the limit that the instrumental noise time stream, n, is stationary, Gaussian, and is well-characterized by a white frequency 
spectrum, the optimal map is obtained by summing all time samples di and decorrelation matrix elements falling in a pixel p. 
Assuming that the scan strategy, instrument, or channel combination provides modulation of the angle ij/, the matrix is nonsingular 
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Fig. 7. The power spectrum of sum and difference time streams from a typical 145 GHz PSB pair observing (apparently) unpolarized atmo- 
spheric fluctuations during the austral summer in Antarctica. When measuring a small polarized signal buried in a large unpolarized background, 
the high degree of common mode rejection of the PSBs makes them naturally suited to an analysis of the sum and difference time streams, as 
described in Section l331 



and the best estimates for /, Q, and U are then obtained by inverting the coadded (3x3) decorrelation matrix at each pixel. This is 
the polarized analog to a naively coadded temperature map. 

The situation becomes markedly more difficult in the presence of noise with nontrivial statistics. The solution that is optimal in 
the least squares sense is again given by Equation[^ with the understanding that now the data consist of the linear combinations 
defined by Equation|2] We now turn to the problem of finding the solution to Equation[^for polarized data, in the presence of 
noise with unknown statistical properties, using channels of varying sensitivity and polarization efficiency. In this general case, 
the data are treated in the following way: estimates of the left hand side of the noise only version Equation|^and the Stokes 
decorrelation matrix are generated for each pixel. 



ni Jj Ci 

y "i Jj Si ) 



(22) 



where the n,- are the elements of signal subtracted time stream. Likewise, for the decorrelation matrix one calculates 



( 1 rjC; TjSi 

2 2 2 

- yf i jjCiSi 

1 1 

- - lYi , 



(23) 



One then obtains an estimate of the corrections to the Stokes parameters /, Q, and JJ maps for an iteration k, by inverting M 
at each pixel. 



(24) 



allowing one to iteratively obtain a solution for the maximum likelihood maps of each Stokes parameter. 



3.5. Sum and difference time streams 

An alternate approach to signal and noise estimation involves operations on the sum and difference of the calibrated time streams 
from bolometers within a PSB pair. This has the numerical advantage of isolating the temperature and polarization terms in the 
numerical inversion of Equation|23] This approach takes full advantage of the high degree of common mode rejection of the PSB 
design, which is illustrated in Figure^ The advantages of this approach, which will be discussed in more detail in Section l33l 
are obtained at the cost of suboptimal noise weighting of the channels within a pair. 

We may represent a sample, /, of a single detector as the linear combination of the sort, 

= / + yi (g COS 2tf/i + U sin 2(A,) . (25) 



Assuming that the channels are properly calibrated, the sum and difference of the signals from a PSB pair may be written as, 

^Si = ^(si + S2)i = / + ^CaiQ + ^(iiU) (26) 
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= - S2)i = ^CaiQ + -piU) , (27) 
where we have defined the angular coefficients 

"^a/ = yi cos 2(/'i,- ± 72 cos 2i/'2i (28) 
*j6/ = 71 sin2(/ri; + y2sin2i/^2/ (29) 

in terms of the independent variables, ij/ifj, where A: = { 1 , 2) identifies the channel. Recall that for a PSB pair the angular separation 
of the channels is A ^ 90 + 2°, however this treatment in no way requires that to be the case. 

Following the prescription of Section l33l one generates linear combinations of the differenced data, 



srai\_ 1 / -aj -arf3i\( Q 



srPi] 2\-ar/3i -0] j\U 



(30) 



As before, one builds up information about the Q, U decorrelation matrix through the combination of channel pairs, as well as 
modulation of the angular coverage, if/. In this regard we have 

A' ■ 

' ' pati s / _ _ \ 

where the time-streams n, represent the polarization subtracted difference data. The 2x2 decorrelation matrix is, therefore, 

and we note that we are now using suboptimal weighting of the pairs to generate corrections to the polarization map. Note that, 
for A ^ 90°, the quantities "a and have opposite parity, so that when averaged over a large sampling of ifr, the off-diagonals 
of Equation |32] are small. Once the corrections to Q and U are obtained, one may substitute them in the sum for ^5 to solve self 
consistently for /. 



3.6. Polarized cross-linking 

The iterative map-making methods described in Sections l3 .41 and l331 result in a self-consistent estimate of the signal and noise 
from the data that is "optimal" in the least-squares sense. However, instrumental effects can introduce correlations in the signal 
contribution to the time domain data that limit the fidelity of the recovered Stokes parameter maps. ^ 

In simulations using noise correlations to process the signal only time streams according to Equation^] these effects appear 
as a residual between the input and recovered Stokes parameter maps. The spatial morphology and amplitude of these residuals 
depend on the amount of cross-linking in the scan strategy, the degree of polarization modulation, as well as the method used 
to decorrelate the /, Q, and U parameters from the time stream. While these residuals do not introduce a bias to the pseudo-Cf 
estimates of the power spectra, they do contribute to the signal covariance of the map, and therefore degrade the sensitivity of the 
Monte Carlo approach relative to optimal methods.^ 

The fidelity of the recovered Stokes parameter maps is an important consideration for the design of scanning polarimeters; 
the statistical depth of the survey determines the level at which these instrumental artifacts must be controlled. To investigate 
these effects, we generate signal-only simulations based on the B03 observation strategy and the measured B03 noise power 
spectra. The B03 cmb data consist of a deep region and a shallow region, representing the extreme cases of possible observation 
strategies avaliable to Antarctic LDB payloads. The B03 scan crosses each pixel in the deep survey over may timescales and at 
many different orientations (due to sky rotation), while the pixels in the shallow survey are not well sampled. 

Using the Healpix synfast facility 1 12|, we generate a noise-free polarized cmb sky, pixelized at 3.4' resolution, from a concor- 
dance ACDM model. We then simulate three polarization modulation schemes to compare with the nominal B03 modulation (i.e. 
sky rotation alone). Each time-domain simulation includes the nominal sky rotation in addition to that which would be achieved 
with a rotating half-wave plate. We model the follo wing modulation schemes, which are representative of those proposed by 
balloon borne and terrestrial bolometric polarimeters 114 ill . 

1. 22.5° steps of the polarization angle each hour. 

' Common examples of such instrumental effects include the impact of the AC coupling of detector outputs, variations in the noise spectra be- 
tween detectors, scan synchronous noise, and limited accuracy of estimates of the low frequency noise. These noise estimates are fundamentally 
sample variance limited by the finite period over which the noise can be considered to be stationary. 

* The residuals contribute directly to the effective transfer function for the temperature and polarization spectra (the T'f discussed in Hivon 
et al. and Contaldi et al. (lo|l'). 
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Fig. 8. A signal-only simulation, showing the residuals between the observed and input polarization (in this case, Stokes Q). The differencing 
method of Section 1331 f /eft panel) is more robust to common mode effects than is the more general method of Section IT^ ( right panel), 
especially in regions where the crosslinking is poor. This is due primarily to the correlations that are introduced by the preconditioning of the 
TODs (essentially a highpass filter at 20 mHz) and the features in the noise kernel A'"', which introduce path dependencies to the observed /, 
Q, and U parameters. 



2. 22.5° steps of the polarization angle at the end of every scan. 

3. Continuous rotation of the polarization angle of each PSB at 350 mHz.' 

We observe the simulated sky with each of these polarization modulation schemes and create a noise-free time ordered data 
set, s, for each. We then solve for the signal part of the general least squares map using the B03 inverse noise filters, A^"', 
according to Eauationll6l 

in = (a^A^^'a)"' A^A^^'s 

The noise kernels, ', are smoothly truncated below 70 mHz.'" We decorrelate the Stokes /, Q, and U parameters at 6.8' 
resolution, using both the general (3x3) method and the PSB sum/difference (2x2) method, and compare the resulting polarization 
maps with the input sky. In the case of the former method, the residuals in the Q and U maps contain contributions from the finite 
resolution of the pixelization as well as the correlations introduced in the course of making maps from the time-ordered data. The 
difference time streams do not contain the relatively large unpolarized contribution, and therefore are far less susceptible to these 
pixelization effects. 

Before considering the effects of the polarization modulation, we first investigate the benefits of exploiting the common 
mode rejection of the PSB pairs through the analysis of the difference time streams. We compare the residuals resulting from the 
application of the general method of Section l34l and from that of Section l331 In Figure|8lwe show the qualitative improvement 
in the fidelity of the reconstruction that results from the analysis of the difference time streams. It should be noted that, even 
for the general method, the sky rotation of the nominal B03 scan provides a degree of modulation that is sufficient to reduce the 
residuals to a level well below that of the instrumental noise in the B03 maps 0(35jl; the use of a waveplate in B03 would 
not have significantly improved the accuracy of the polarimetry. Furthermore, the direct difference method of Section l331 is less 
sensitive to the limited cross-linking of the nominal scan than is the general polarization decorrelation method of Section l3r4l 

While the design of the PSBs is naturally suited to the sum/difference approach, scanning experiments sensitive to a single 
polarization (such as Ebex |41] and Spider |39]) are not able to exploit the common mode rejection that is intrinsic to the design 
of the PSBs. A scheme for polarization modulation is therefore a highly desirable feature in singly polarized systems. 

In order to illustrate the effect of the polarized cross-linking on the fidelity of the reconstructed signal, we show in Figure|5] 
the residuals that result, for a particular cmb realization, from the general (3 x 3) approach of Section IT^ for the nominal B03 
scan, and for each of the three modulation schemes listed above. The fidelity of the signal reconstruction improves in proportion 

' We choose this modulation rate to be as fast as possible, given the Boomerang scan and sample rates. 
The lowest frequency that can be reliably recovered clearly affects the amplitude of the residuals, especially for the temperature fluctuations 
which have significant power on large scales. However, the relative benefits of the scanning strategies outlined above are generally insensitive 
to the exact value of the minimum frequency. 
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to the rate of modulation. The improved cross-linking randomizes the path dependencies of the observed Stokes parameters that 
are introduced by the processing of the time ordered data, namely Equationll6l 

As shown quantitatively in Figure^] the largest residuals in the Q and U maps (which, in the B03 example, occur on rel- 
atively large scales) can be significantly reduced by a modest degree of polarization modulation. Nevertheless, for the effects 
that have been included in this simulation, generating maps using the sum and difference time streams of the PSBs is nearly as 
effective at minimizing the residuals in the Q and U maps as the use of a half wave plate. It is important to note that beam asym- 
metries, instrumental polarization, pointing errors, and calibration uncertainties are examples of effects that are not included in 
these simulations. Each of these effects are mitigated by the use of an idealized modulation scheme, but not by the sum/difference 
method of Section l331 

4. Summary and Conclusions 

We have described in detail the design and performance of the Polarization Sensitive Bolometers (PSBs) which have enabled 
the first generation of successful bolometric cmb polarimeters. This discussion outlines the instrument parameters which must be 
characterized to accurately decorrelate the Stokes /, Q and U parameters from the time ordered data of a PSB. The design of 
the PSBs provides a high degree of common mode rejection that can be exploited in the analysis to minimize susceptibility to 
various instrumental effects that can potentially limit the fidelity of the recovered polarization. Simulations of PSB data including 
realistic instrumental effects illustrate the benefits of analyzing difference time streams, as well as various practical schemes for 
modulating the polarization signal. 
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Fig. 9. The Joy of Crosslinking. The residual signal in the Stokes Q ( left column ) and U ( right column ) parameter maps generated 
using the general algorithm of Section for increasing levels of polarization modulation. At top, the only modulation is that 
provided by sky rotation. The middle two sets are the residuals obtained when stepping the half waveplate by 22.5° (Q — > U) each 
hour and at the end of each azimuth scan, respectively. The bottom row shows the fidelity achieved with a waveplate spinning 
continuously at 350 mHz. The remaining residuals are dominated by pixelization effects. 
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Fig. 10. Histograms of the Stokes QjU residuals, using 6.8' pixels, for various modulation schemes. Figure|9|shows that the errors are largest 
on large scales. The residuals in the last row are dominated by pixelization effects. 
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Appendix A: Noise in bolometric receivers 

General noise properties 

Noise in bolometric receivers originates from several independent sources, including contributions from the readout electronics, 
the detector, and the intrinsic fluctuations in the optical background power These noise sources are independent processes and 
their contributions add in quadrature to the total noise of the system. 

The contribution of each of these components to the time ordered data are filtered by the transfer function of the one or both 
of the detector and the readout electronics". The voltage noise (that is, the Johnson, JFET/amplifier noise, and the product of 
current noise with the series impedance) of the system, n,., is filtered only by the transfer function of the readout, Z(/). The 
photon noise, «y, and phonon noise, rtc, are filtered not only by the bolometer voltage responsivity, S (/), but also by the readout. 
The bolometer transfer function, S , is typically that of a single low pass filter, or a cascade of two such filters. The details of the 
readout transfer function, Z, vary, but always includes an anti-aliasing filter which strongly attenuates frequencies well below the 
Nyquist frequency of the analog-to-digital converter (ADC). The raw data are composed of the signal from the sky, s, and the 
various noise contributions, convolved with the bolometer and readout transfer functions, 

d-Z® [5 ® {s + ny^nc^ + n,,j . (A.l) 

One of the first stages of analysis involves the deconvolution, and then the de-glitching, of this raw detector time stream. The 
deconvolution is normally accomplished in the Fourier domain by dividing the product of the bolometer and electronics transfer 
functions, Z' - Z S , from the raw data, d. While this results in a time stream that is characterized by a signal component with 
a uniform calibration in the frequency domain, the contribution of the voltage noise is biased according to the detector transfer 
function, — > HylS . Because the bolometer transfer function has the form of a low -pass filter, the deconvolved time stream 
generally exhibits "/-noise" in proportion to the time constant of the bolometer and the amplitude of the voltage noise component. 

Figures and Is] are examples of the power spectrum of a time stream prior to deconvolving the system transfer function. 
Figure |6l shows the power spectrum of the deconvolved data. A precise knowledge of the system transfer function is necessary 
to avoid the introduction of instrumental artefacts in the recovered signal. In the case of CMB studies, such an error generally 
results in a bias in power spectra derived from the maps iEtIi . 

Contemporary bolometric receivers, even those operating in the low background environment provided by balloon and orbital 
payloads, are designed to achieve background limited sensitivities. In these receivers photon noise is a major, if not dominant, 
contribution to the total noise in the system. In the following section we examine aspects of this photon noise contribution, Uy, 
including a derivation of the (low level) noise correlations expected between detectors in a PSB pair resulting from fundamental 
properties of fluctuations in the thermal background radiation. 



Photon noise 

A fundamental limitation to the sensitivity of any receiver (band-gap, coherent, or bolometric) derives from the intrinsic tem- 
poral fluctuations in the optical, often thermal, background radiation. The noise properties of thermal background radiation, or 
photon noise, differ greatly between radio, sub-millimeter, infrared, and optical instrumentation due to their vastly different op- 
erational regimes of photon occupation number. Photons satisfy Bose-Einstein statistics, and therefore the occupation of a mode 
of frequency v is 

for a thermal background of temperature, T . 

The ratio kjh - 20.8 [GHz/K] sets, for a given background temperature, the frequency for which average occupancy is above 
or below unity. At radio wavelengths astronomical instruments typically enjoy background levels of order 10 K, with the min- 
imum background limited by the cmb monopole at 2.728 K. At higher frequencies, atmospheric loading and thermal emission 
from the instrument tend to dominate the background, and are typically ~ 30 - 100 K for terrestrial telescopes. Therefore, instru- 
ments operating at frequencies above ~ 100 GHz have occupation numbers of order unity, while receivers at lower frequencies 
tend to have very large occupation numbers, n ^ kT/hv. In the low n regime, photons can be thought of as arriving at the detector 
sporadically. The photon noise in high frequency (> 100 GHz) instruments with low backgrounds can therefore be expected to 
largely satisfy Poisson statistics, where one expects fluctuations on the mean to scale roughly as V^V. 

" For contemporary bolometric instruments like Planck HFl, the readout electronics include a cold JFET amplifier, ambient temperature 
amplifier/bandpass filters, and an anti-aliasing/data acquisition system. In practice, the contribution to the noise of everything except the JFET 
amplifiers and low noise preamplifier are neghgible. 

'■^ In some instances, the electronics are AC coupled, meaning that the low frequencies are strongly attenuated. The advent of low-cost, 
high-resolution ADCs has made this feature less common. 



20 



W. C. Jones et al.: Instrumental and Analytic Methods for Bolometric Polarimetry 



Hanbury Brown and Twiss were the first to complete a rigorous analysis of noise correlations in photons 
The topic has been continually revisited in the fifty years since the first published work, and is still relatively un-advertised among 
many instrumentalists and observers alike. Therefore, we go through the analysis in detail. 

Following Zmuidzinas |53J, we can write the covariance matrix describing detector outputs in all generality 



4 



l^lj dvBij{Bji + dij) , (A.3) 
where we define the power coupling matrix 

Bij = hvJ^S ikS + Cij , (A.4) 

k 

and the internal noise term, 

I hv e" + 1 

C,v = (I-SS%y^;^. (A.5) 

Here, S ij is the standard scattering matrix, which couples an output amplitude to the inputs at each port of a network, as in Figure 
lA.ll It is defined by a, - Yi±§ U^i- Eauation lA.5l x = hv/kTs, where Ts is the thermodynamic temperature of the system 
S , and the in Equation IA.4I are the occupation numbers of the modes at port k. We take ports and 1 to label the two input 
polarization states, and let ports 2 and 3 label the two bolometers in a PSB pair. For simplicity, we assume that n2 = "3 = (i.e., 
the detectors are extremely cold with respect to the background), and that the input populations - n\ - n{v, Tioao) imply no 
net polarization in the background. 

The internal noise term, C/j arises as a result of losses in the system. Any mechanism causing loss implies a thermal noise 
contribution, c,, to the outgoing signals 

fl; - 'LjS ijbj + Ci 

which depends on the temperature of the lossy component. 

In an ideal lossless network, the system's thermal noise term will vanish since S is unitary, (I - S 5 ' = 0. In this case, 
the only nonzero terms in the scattering matrix are 5 20 - 5 31 - 1. Since the only nonzero populations are no,i = n, the 
covariance matrix contains only terms with B20 = ^31 = « ■ hv- Under this assumption the detectors' noise is uncorrected, and 
the autocorrelations satisfy 

aji = J dvnin + 1). (A.6) 

Therefore the Icr uncertainty in the incident power due to intrinsic background fluctuations is 

hv / Av I 

■ ^jr]n(J]n + 1) . (A.7) 



^photon 



T] \ T 



Note that in the above we have explicitly included the optical efficiency, rj = |52ol^ = |53ip. This result differs from the familiar 
Dicke radiometer equation describing coherent receivers, 

hv /Av 

O-photon = \ (W + 1). (A.8) 

77 V T 

While Eauation lA.Sl is the limiting form of Eauation lA.7l for large occupation number «, it is instructive to derive Eauation lA.81 
from Eauation lA.3l 



Example 1 : The Dicke radiometer equation 



The scattering matrix for an idealized coherent receiver with perfect isolation contains a single nonzero term, 15 loP = G, where 
G is the gain of the system. An amplifier can be thought of as a population characterized by an inverted distribution of energy 
levels, such as that found in a maser or laser Such systems are conveniently described in terms of a negative temperature. As 
T — > 0, the sign of the Cij from Eauation lA.Sl is reversed. The only nonzero element of C is Cn = G - 1, and therefore 
Bii = Gn + G - 1. Application of Eauation lA. 31 gives 



^ J dyG'[n + l-iGr'][n + l] 



Av 



(hvf 



in + If 



(n+l)l 



(A.9) 
(A. 10) 
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Fig. A.l. The scattering matrix of a four port network. A polarization sensitive bolometer can be modeled by such a network. 



In the limit that G is significantly larger than unity, the second term becomes negligible. Referencing the noise to the input, we 
recover the (lossless) Dicke radiometer equation, Eauation lA.81 

(Til = "I^a/ — {n + 1). 
\ T 



Example 2: Polarization sensitive bolometers 

A dual polarized, single-moded receiver (coherent or bolometric) is completely described by a four port network. Polarization 
sensitive bolometers and coherent receivers using orthogonal mode transducers (OMTs) are two examples of such systems. We 
now derive the photon noise properties of a PSB pair 

The action of the network, S , is that of an imperfect polarized beam splitter, with two inputs and two detectors. The PSBs 
(two of the four ports, labeled say, as numbers 2 and 3) are assumed to be at cryogenic temperatures, and therefore contribute 
negligibly to the photon occupation number. Therefore, the entries in the scattering matrix relevant to the observed photon noise 
are limited to the lower left quadrant, namely 520 = y, 5 31 - y\S2\ -5, and 5 30 = S' . Here the parameters y and 5 describe the 
efficiency of transmission of the copolar amplitude and the crosspolar amplitude, respectively, and in practice j » 6. 

Only the lower right quadrant of SS ' is nonzero. 



{S s 022 = + S6' (S s 1)23 = 5 (r + r') 
(55 ^)32 = 5' (r + r') (SS'^hi = r" + S6' 

The power coupling terms, Bij, of interest are given by 

B22 = (|52ol'no + |52il'ni) hv + C22 

- (r^ "0 + 6^ rii + [I - (j^ + 6d'y\ «£.) hv 
Bn - (S20S30 no + S21S31 ni) hv + C22 

— {yd' riQ + y'6 rii - 6 (y + y') tic) hv 

where we have written the thermal contribution of the network 



(A.ll) 



(A. 12) 

(A. 13) 
(A. 14) 
(A. 15) 



1 1 

tlr = . 

2 - 1 

In the case of PSBs, the source of the modal coupling are the detectors themselves and/or the optics. In the case of the detectors, 
they are extremely cold compared to the background. For Boomerang, the reimaging optics and filters are also cooled, and have 
low emissivity. Therefore, we assume the thermal noise contribution of the network, ric, is very small compared to the background 
populations, n,. Furthermore, we assume that the background is isotropic, i.e., no = "i = «■ The covariance of the photon noise is 
then fully described by 



(tI = (hvf^ [ (/ + + (r^ + 6^)n\ 



ajj = {hvf—[(2y5fr?\. 



(A. 16) 
(A.17) 
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The autocorrelation, Eauation lA.161 contains terms proportional to both rp- and n. The former is commonly referred to as the Bose 
contribution, or as a "photon bunching" term. Equation lA. 1 7l shows that correlations between devices are proportional only to the 
Bose contribution, implying that PSBs operating under higher background loading conditions will exhibit a higher proportion of 
correlated noise than the same instrument operating in a lower background. For an idealized system, in which the polarization 
leakage S is zero, the covariance between detectors vanishes since the two linear polarization states are statistically independent 
of one another 

In practice, we estimate the total optical background power, Q, arising from the cmb, atmosphere, the telescope, and emission 
from within the cryostat. For simplicity, this optical background is treated as having originated from a single thermal source at an 
effective temperature Tju - Q/TjksAv. The noise equivalent power from the background fluctuations is then given by Equation 
|A36l 

NEP2,„,„„ - 2hy e (1 + 77 «(r«y)) . (A. 18) 

This is, of course, only approximate as we do not treat the background sources independently. It is often the case, however, that 
a single thermal source contributes the majority of the background optical power 



Appendix B: The Jacobi method 

We outline the application of the Jacobi method to the problem of mapmaking from scanning experiments, much of which can be 
generalized to other iterative algorithms such as the method of preconditioned conjugate gradients. We loosely follow the more 
complete discussions of the topic which can be found, for example, in Acton yj, Barrett et al. |^, Young 115 ill , and Press et al. 

The Jacobi method is a robust numerical method of solving a set of linear equations. Ax = b, for which the matrix A is (or 
can be arranged to be) diagonally dominant. The great strength of the Jacobi method is that, subject to this requirement, it is 
guaranteed to converge although it may do so relatively slowly. Given a trial solution, one may estimate a new solution without 
inverting the matrix A simply by solving for each component, x*^', given an estimate of the values {x*). 



It is often convenient, and advantageous from a numerical point of view, to write the above in terms of a correction to the 
previous iteration. 



X?+' = X* + 5xf+' 



where 



5x^^77 A,' 



(B.l) 



Here we have inserted a convergence parameter 77 < 1, which may be tuned to aid the convergence of the algorithm. In the limit 
that A is diagonal, the optimal value isT] - 1 . Generally speaking, the larger the off-diagonal terms become, the lower the optimal 
value of T]. Clearly the diagonals of A must not be near zero. Furthermore, as can be seen from Eauation lB.il the solution will 
diverge if the absolute value of the sum of the ofF-diagonals is greater than the diagonal element of each row. 
As an example, consider the following linear system: 



' 5 


-2 


1 ^ 




'-1 ^ 
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-7 
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X = 





1-2 
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6j 




[ 1 j 



(B.2) 



Setting ?7 = 1, and using the above procedure results in the following sequence of solutions: 

xo = ( 0.000, 0.000, 0.000) 
xi = (-0.200, 0.000, 0.167) 

X5 = (-0.291,-0.187, 0.102) 



Xoo = (-0.300,-0.200, 0.100) 
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This example converges to twelve significant digits after 40 iterations, largely independent of Xq, the trial solution. The rate of 
convergence does not scale strongly with the array size, so the solution is an efficient way of solving large systems of equations. 

A minor modification to the above procedure results in the Gauss-Seidel algorithm, for which the estimate for each value 
x'l^^ incorporates the most recent estimate of the parameters {x*^' )j<,- instead of the set of values from the previous iteration. This 
procedure is less numerically robust, but tends to converge more rapidly than Jacobi iteration. 

The application to Equation^]is clear; the Jacobi method provides a robust method of solving for the general least squares 
map. Equating Equations^]and Eauation lB.ll we find the correspondence, 

A ^ C^' = (A^N^'A) 

X — > in 

b ^ A^N 'd 

Recall that the matrix A, which appears on the right hand side, is the pointing matrix and should not be confused with the general 
linear system described in Eauation lB.il The algorithm we use for calculating the correction to an estimate of the least squares 
map, in*, is simply 

5m*+' oc diag(A^N"'A)"'A^N"'(d-Am*) . 

In practice, one can simultaneously solve for the noise covariance matrix of the data, N. This typically results in slightly slower 
convergence of the algorithm than when using a fixed noise estimate. 



